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^^ Abstract. The vector Burgers equation is extended to include pressure 

gradients and gravity. It is shown that within the framework of the Cole- 

Hopf transformation there are no physical solutions to this problem. This 

^ ' result is important because it clearly demonstrates that any extension of 

^D ' Burgers equation to more interesting physical situations is strongly limited. 

O ' 



1. Introduction 

We solved Burgers equation in three dimensions for arbitrary orthogonal curvilinear coordinate systems 
[1]; hereafter, NSM (also see [2]). The analytical solutions suffered from the well-known deficiencies of 
Burgers equation, which does not account for gravity and pressure gradients, nor does it conserve mass. 
In this paper, we have attempted to extend Burgers equation by including gravity and either isothermal 
or adiabatic pressure gradients. The basic idea of our approach is to split the momentum equation into 
the vector Burgers equation and the equation which contains all remaining terms, and to solve these two 
equations simultaneously. We show that the latter can be done only when a special auxiliary condition is 
satisfied. Because the solutions to the vector Burgers equation are well-known [1], the main problem is to find 
physically acceptable solutions to the auxiliary condition that also satisfy the continuity equation. We have 
attempted to generate steady-state, viscous solutions by constructing them from the solution of Laplace's 
equation in such a way as to satisfy the continuity equation for compressible flow. Our apparent success is 
vitiated when we find that there are no physically acceptable solutions to the 2-D viscous case that can be 
derived from the potential flow problem and produce physical pressure gradients. This is an important result 
because it shows that the obvious extension of the vector Burgers equation to more interesting physical cases 
is strongly limited. Still, the solutions arc of interest in the same sense that the analytic solutions of the 
vector Burgers equation are interesting. 

We begin our presentation by summarizing the assumptions and solutions of the vector Burgers equation 



presented in paper [1] and then describe in detail our new results. 



2. Basic Equations 



The vector Burgers equation can be written in the foUowin form: 



^- + u • Vu = i^V-^u. (1) 

ot 



We proved that the solution of this equation can be expressed as a function of the solution of the diffusion 



equation: 



se ^2. 



-m - ^^^' (2) 



namely, 



"=-§^Q- (3) 



This result is correct for arbitrary curvilinear coordinate systems and is a subset of the more general so- 
lutions to the tensor Burgers equation [2]. We would like to be able to extend our previous solutions to 
include pressure gradients and gravity, as well as satisfying the continuity equation. Then we must solve the 

momentum equation, 

du ^ ^ VP GM o^ - 

-—+u-Vu^ hV \-iyV^u+f (4) 

ot p r 



where / is a general conservative force per unit mass so that 



/ = V*/. (5) 
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and /9, P, G, M are density, pressure, the gravitational constant, and the mass of the object exerting the 

gravitational pull, respectively. The continuity equation is written as 

l^+V-M^O. (6) 

In order to evaluate the vector Laplacian for arbitrary orthogonal coordinate systems, wc must use 

V2u = V(V-m)- Vx (V xu), (7) 



and for the inertial term 



t?- Vu = V— -u X (V X u). (8) 



Then eq. (4) greatly simplifies as long as 

u X (V X m) = vV X (V X u). (9a) 

This appears to be an equation that represents quite general flows, but the only solutions are irrotational. 
This can be seen by taking the divergence of both sides of eq. (9a) and using 

V -{Ax B) ^ B-{V X A)- A-{V X B). {%) 

Therefore we require 

Vxu==:0, (10) 

and 

u = V*. (11a) 



We also use 



p \J p 



(116) 



so that all terms in eq. (4) are the gradient of a scalar. The integral of the pressure term can be done whenever 
P = P{p), but we will solve eq. (4) in the limit of an isothermal pressure gradient with the ideal gas law 



(P = a^p, where a is the isothermal sound speed). Then the pressure gradient term may be written as 



VP 



= a^Vlnp. 



(12) 



Now, using eqs. (9a), (11a), and (12), eq. (4) simplifies to 



9$ (V$)2 
Ih^ 2 



a^ In p - 



GM 



jyV^$-$/ 



(13) 



This implies that the integral of eq. (13) is only a function of time, so that 



5$ (V$)2 
'dt^ 2 



CM 
a^lnp- i/V2$ - ^f = E{t). 



(14) 



The following vector identity is useful and leads to an important form for V^$: 



va^llve, 



(15a) 



which is true in any orthogonal coordinate system. Taking the divergence of both sides with a = $ and 



using eq. (15a) again with a = ^ yields 



v^* = 0,ve)' + *v»e, 



dQ 



(156) 
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which, again, is vahd for any orthogonal curvihncar coordinate system. Substituting eq. (15b) into (14) 



a^ Inp ----$;- W ;^(Ve)^ + — V^O = E{t). (16a 



dt 2 ^ r ■* Vde2' ' de 



Now cq. (2) guarantees that the first and last terms on the left-hand-side will cancel so that 



(ve)2 



1 fd<i>V £^ 



2\dQ) ''de^ 



a'^lnp- — -<!>f = E{t). (166) 



An important aspect of our work in NSM was that the remaining terms in the equation similar to eq. (16b) 
were proportional to (VO)^ so that this factor cancelled out, leading to an integrable equation for ^. It is 
clear that the same procedure will work here provided that the following auxiliary condition is satisfied: 



CA'T 

aMnp--:^ $.-S(t)=±C2(Ve)2. (17a) 

r 



It is easy to show that C must be found from the solution of 



-"<«)^^§-KS)'- 



but this is not particularly helpful. We choose to examine possible solutions for the case of constant C 
partly because C=0 corresponds to the vector Burgers equation solution discussed by NSM, and we are 
attempting to find physical extensions of the Burgers equation solutions. We are also physically motivated 
by noticing that the particular choice of C^ = —1/2 in eq. (17a) reduces the auxiliary condition to the 
Bernoulli equation which must be solved for the pressure distribution implied by a potential flow model (the 
steady limit of the diffusion equation). The difference is that the auxiliary condition is for compressible 



7 
flow. Additionally, as we will also show, the case of constant C generalizes known velocity laws derived from 

the solution of Burgers equation. It is for these reasons that we examine the solutions of the irrotational 

momentum equation subject to the constraint of eq. (17a). 

The differential equation that results from substituting eq. (17a) into eq. (16b) (the analog of eq. (7a) in 
NSM) is 



We solve for ^ and use 



^ ^2^^^^^- ^''^ 



" = V$ = -^Ve, (19a) 

dfc) 



to find 



de 



V2C ta.n ( -^e + Ci 



(196) 



u = %/2Ctan(^e + Ci j Ve, (19c) 

where Ci is a constant of integration. This solution corresponds to the minus sign in eq. (18) while 

w = %/2Ctanh(^e + Ci j Ve, {19d) 

corresponds to the positive sign. The first solution seems nonphysical because of the singularities, although 
solutions could be generated by restricting the argument between and 7r/4. But, as we will show shortly, 
the existence of any physical solutions require choosing the hyperbolic tangent solution. It is also known 
[3,5] that the limit of eq. (19d) with V8=constant is the steady solution to the 1-D cartesian Burgers 



equation for an initial step function profile with a sharp but finite diseontinuity in u at x=a whieh rounds 
off asymptotically in time to a steady-state shock layer profile of thickness v /uq. 



u = — «o tanh 



Mo (a; — a) 



2v 



(19e) 



The tangent form is also known to be a steady-state solution of the cartesian Burgers equation [5] . 

The obtained solution, eq. (19d), leads to a mechanical energy equation which is found by taking the dot 
product of eq. (19d) with u and then eliminating (V0)^ in eq. (17a). This gives 

,,2 ( C \ CM 

— coth^ f-^e + Ci) + a^ In p-<Pf^E{t). (20) 

Note that although this equation is similar to Bernoulli's equation, eq. (20) is for nonsteady, viscous flow. 
The hyperbolic tangent solution can be seen as the appropriate form because the internal energy/gram 
will decrease if the speed increases. The form of the mechanical energy equation is strongly reminiscent of 
nonlinear decay models [1]. A singular source of mass at the origin requires solutions for cothO that begin 
at finite radius close to the source, but then the motion decays so that the kinetic energy/mass approaches 
zero as t^ oo when E(t)^ 0. 

We must now impose mass conservation on this model, and, to this end, we find a useful form for the 
continuity equation. Using eq. (11a) and further eliminating V • u by using eqs. (2) and (15b) yields: 



^+.-V lnp+— (Ve)^ + --^0. (21) 



To evaluate the continuity equation, we use streamline coordinates 



e. = ^ (22) 



3. Steady Solutions 

We further limit physical space to steady solutions so that cq. (2) reduces to Laplace's equation. Then eq. 
(19d) requires that VO is zero perpendicular to streamlines and V0 is ^. The problem reduces to finding 
that subset of the solution of Laplace's equation that satisfies the auxiliary condition and conserves mass. 

Identifying V0 with a velocity derived from a potential flow solution 



V • Upot = 0. (23a) 



V • {pu) ^ 0. (236) 



so that 



pu == PpotUpot + V X A. (23c) 



The following details show that V x A is zero. 



The steady continuity equation is derived from eq. (21) 



d , d^<^ fdQ^^ 



"a^^^^+^U) -'■ ^''^ 



Dividing by ^ and using 
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(for solutions with non-zero velocities), eq. (24) integrates exactly with respect to Q to 



,% = '.. i^o) 



where ki might be a function of spatial coordinates. Substituting eq. (26) back into (19a) 



pu = kiVe, (27a) 

or, if the original solution of Laplace's equation is formulated in terms of the velocity, 

pu = kiUpot, (276) 

and fci is seen to be the original incompressible density for the potential flow solution, Ppot, i-e., a strict 
constant. An important aspect of this method is that the new steady, compressible solutions are created 
from the old potential flow (i.e., incompressible) solutions, provided that the boundary condition on the 
compressible solution is formulated in terms of the mass flux. 

Substituting the potential flow speed back into eq. (20) yields another form for the steady mechanical energy 
equation: 



aHnp~ — ~'^f + C^ul^t^Eo. (28a) 

Setting <!>/ to zero and choosing 

C - I, (286) 
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we get 

a^lnp- — + % = i?o. (29) 

r 2 

The auxiliary condition is now seen to be the Bernoulli equation that must be used to solve for the pressure 
distribution implied by the potential velocity. However, we use eq. (19d) to solve for the next higher-order 
solution for the velocity related to this variable pressure distribution, but including viscosity. 

With this value for C, we obtain 

u = - tanh ( ^ + C*! ) ^pot, (30a) 

and 

p = Ppot coth ( ^ + C*! ) ' (30fc) 

The last part of this solution is to solve for the subset of physical potential flow solutions. Setting gravity 
to zero and substituting eq. (30b) into (28a), we find 



ul^^ = 2Eq - 2a^ In coth ( -^ + Ci ) . (31a) 



Replacing the J dP/ p with the adiabatic term, (7 — \)P/^p, and substituting 7 — 5/3, and then using eq. 
(26) gives the related equation for the adiabatic case: 



■2 = 2Eo - Sfecoth^/^ ( ^ + Ci ) ' (31^) 



where we have also used the adiabatic equation of state. 



p = bp'^. (31c) 
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The potential flow solutions are limited to those whose speed is constrained to obey eqs. (31a, b). 

We use complex variables to show that there are no physical solutions that simultaneously satisfy both 
Laplace's equation as well as the constraint on the speed, eqs. (31a, b). Writing the complex potential, F(z), 
as 

F{z) = Q{x,y)+i-^{x,y), (32a) 

with 



z ^ X + iy, (326) 



then 



dF_ 
dz 



W{z) ^ — ^v,- iVy, (32c) 



where ^ is the stream function, and W(z) is the complex velocity. The utility of this approach is obvious 
because 

WW*^vl+vl, (33) 

where W* is the complex conjugate of W, which allows the velocity components to be read from eqs. (31a) 
or (31b). The method of complex variables requires that these velocity components must be derived from 
an analytic potential function. In the following, we demonstrate that this is not possible. 
Consider the general form of eqs. (31a,b): 



(ve)' = /(e). (35) 
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This equation implies that the gradient of (VO) must be parallel to VQ, or that the Jacobian of 8 and 

(VO) must be zero. This is a necessary condition for solutions to eqs. (31a, b) to exist. The Jacobian 

simplifies (in cartesian coordinates) to 



2 2^dv^ o ^'" 



{vl^vl)^-2v,Vy-^^0, (36) 



where we have also used the Cauchy-Riemann equations, 



9$ 9* 



But eq. (36) is the same as 



r*2..dW{z] 



Im[ W*\z) ^ ^ ) =0. (38a) 



One solution is the trivial case of zero derivative. The other is found by defining G(z) to be the reciprocal 
of W(z), so that eq. (38a) may be written as: 



-(ii^)-- w 



Because the denominator is real, we require, except at poles. 



-(^1-- <-) 



But eqs. (37) then require 

dG{z) 
dz 



= 0, (39) 
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so that 

W{z) = -, (40) 

z 

where c is a complex constant. This gives the solution 

8 = ailnr + a2, (41a) 

together with the trivial solution 

6 = aaa; + a4y + as, (41^) 

where the a, are constants. As cqs. (31a) and (31b) are not of this form, there are, then, no two-dimensional 
solutions of the potential flow problem that satisfy the physical constraint on the speed. As these are 
required to solve for physical viscous, compressible solutions through eqs. (29) and (30a), there are no 
physical solutions of this model. This is an important result because it shows that the method of solving 
Burgers equation based on the the Cole-Hopf transformation cannot be used to solve a general (and physically 
more interesting) problem of a flow with gravity and pressure gradients. To sum up, substituting eq. (26) 
into eq. (16b) clearly indicates that (V0)^ is a function of O. But, we have shown that the only steady 
solutions for V8 arc those generated from eqs. (41a, b). Our Colc-Hopf extension of the solution of Burgers 
equation requires that F{Q) must be a function of coth(0/2i^ + Ci), so that there are no physical solutions. 
4. Conclusions 

We have extended the solutions to the vector Burgers equation [1] by a method that splits the momentum 
equation into two pieces, each of which is satisfied simultaneously. This allows us to solve a subset of 
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irrotational flow problems with pressure gradients and gravity. The steady hmit, when combined with mass 

conservation, reduces to specifying a subset of potential flow solutions. These incompressible flows, in turn, 

allow the generation of viscous, compressible solutions that are related to the solutions of the vector Burgers 

equation [1] . We then proved that the subset of two-dimensional potential flow solutions does not correspond 

to the requirements for physical flow derived from real pressure gradients. In this sense, the pressure gradients 

are derived from pseudo-densities. It seems unlikely that physical models could be generated by this method 

using either three spatial dimensions or time-dependent solutions. Still, these results are of interest in the 

same sense that the analytic solutions of the vector Burgers equation are interesting. 
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